Phase Dependent Spectral Variability in 4U1907+09 
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ABSTRACT 

We report on ASCA, RXTE, and archival observations of the high mass X-ray binary 
pulsar 4U1907+09. Spectral measurements of the absorption and flux were made at 
all phases of the X-ray pulsar orbit, including the first spectral measurements of an 
extended period of low flux during two of the ASCA observations. We find that a 
simple spherical wind model can flt the time averaged light curve as measured by the 
RXTE ASM, but does not flt the observed changes in the absorption column or account 
for the existence of the phase-locked secondary fiare. An additional model component 
consisting of a trailing stream can account for the variations in column depth. However, 
these models favor a high inclination angle for the system, suggesting a companion mass 
more consistent with an identification as a Be-star. In this case an equatorially enhanced 
wind and inclined neutron star orbit may be a more appropriate interpretation of the 
data. I 

Subject headings: X-rays: stars — pulsars:individual (4U1907+09) — stars: emission-line, Be 
— stars: mass loss 



1. Introduction 

4U1907+09 was noted as a variable X- 
ray source during the early sky surveys with 
Uhuru (Giacconi et al. 1971), Ariel V (Se- 
ward et al. 1976), OSO-7 (Markert et al. 
1978), and HEAO-1 (Schwartz et al. 1980). 
The latter allowed Schwartz et al. (1980) to 
identify the source with a heavily reddened 
rriy ~ 16.4 star having a broad Ha emission 
line, which they assumed was an OB super- 
giant. A large X-ray outburst observed by 
Ariel V during January 1980 indicated regular 
flux modulation on a time scale of days. Fur- 
ther analysis of all the Ariel V data revealed 
a period of 8.4 days (Marshall and Rickets 
1980). The folded period profile shows a large 
primary peak as well as a smaller secondary 
maximum ~half a cycle later. Long term ob- 
servations by Vela 5B seemed to indicate an 
additional quasi-periodicity at P ~ 41.6 days, 
which was proposed to be due to precession 
similar to that of Her X-1 (Priedhorsky and 
Terrell 1984). 

In 1983, Tenma observed 4U1907+09 for 
an entire 8 day cycle, which resulted in the 
discovery of double peaked 437.5 s pulsations 
(Makishima et al. 1984). Analysis of the 
pulse arrival times allowed an estimate of the 
orbital parameters to be made, which con- 
firmed the 8.38 d period as the orbital period. 
A moderately eccentric (e ~ 0.22) orbit best 
fit the data, although a circular orbit also pro- 
vided an adequate fit. Makishima et al. also 
determined a mass function of ~ ^Mq, con- 
firming the massive nature of the companion. 
Further observations by EXOSAT (Cook and 
Page 1987) implied a mean spin-down rate 
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of ~ 0.23 s yr~ , a variable column density 
(1.5 — 5.7 X 10^^ cm~^), and slightly refined 
the Tenma orbital elements. The observed 
spectra were consistent with a hard power- 
law (a 1.2) with variable low-energy cut- 
off. Balloon observations of 4111907+09 in 
the 20-100 keV range and analysis of Einstein 
archival data are consistent with the EXOSAT 
spectral observations (Chitnis et al. 1993). 

Ginga spectral studies suggested a cy- 
clotron absorption feature at ~ 19 keV (Mi- 
hara 1995; Makishima et al. 1999), which 
has been confirmed by a recent SAX observa- 
tion (Cusumano et al. 1998), along with the 
discovery of the second harmonic at 39 keV, 
suggesting a surface magnetic field strength 
of ~ 2 X 10^^ Gauss. An iron emission line 
at ~ 6.4 keV was also detected, with an 
equivalent width of ~ 60 eV. A pulse pe- 
riod of ~ 439.5 lb 0.6 s is clearly seen in the 
Ginga data during most of the observations 
where the source had an average intensity of 
15-20 mCrab. However 4U1907+09 also ex- 
hibited extended periods of weaker intensity 
(~3 mCrab) accompanied by a reduced pulsed 
fraction (Makishima et al. 1999). 

lye (1986) performed the first detailed 
spectroscopic studies of the companion star to 
4U1907-I-09. He concluded, in contrast to the 
initial assumptions by Schwartz et al. (1980) 
that the spectral type of the companion was 
B2 HI - Ve, and that the two X-ray flares seen 
in some orbits were due to the intersection, 
and subsequent accretion of material in an 
extended disk around the companion by the 
orbiting neutron star. 

van Kerkwijk, van Oijen, and van den 
Heuvel (1989) next undertook extensive spec- 
troscopic studies to determine whether the 
companion to 4U1907+09 is an OB super- 
giant or a Be-star. They considered X-ray, 
optical and orbital data and suggested most 
of the measurements support the supergiant 
nature of the companion, including redden- 
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ing and distance, FWHM, and rotation 
velocity of the companion. In addition, the 
location of 4U1907+09 in the spin-period ver- 
sus orbital-period diagram (Corbet 1984) is 
inconsistent with the trend observed for most 
Be/X-ray pulsars, suggesting instead a wind- 
fed supergiant system. Since it is difficult to 
produce a phase-locked secondary flare from 
the wind of a supergiant star, van Kerkwijk, 
van Oijen, and van den Hcuvel (1989) sug- 
gested the secondary flare may be spurious in 
nature. There is one clear case of a wind- 
fed supcrgiant/X-ray pulsar system, GX 301- 
2, which has a tendency to flare twice per or- 
bit for extended periods of time (Pravdo et 
al. 1995), but the cause of this is not well 
understood. 

Recent observations of 4U 190 7+ 09 with 
RXTE showed highly variable flux with fre- 
quent drops to below the detection threshold 
(in't Zand, Strohmayer and Baykal 1997). 
There is no evidence for increased absorption 
during the drops in their data, suggesting the 
low flux levels are due to a low accretion rate. 
A possible flare that is consistent in phase 
with previously observed secondary flares was 
also seen in the RXTE data. The 2-30 keV 
RXTE spectrum is well fit by a cutoff power- 
law having a photon index a ~ 1.0 — 1.3 with 
a cutoff energy of 13.6 keV, and Nh varying 
from 1— 9x 10^^ cm~^, the lowest values occur- 
ring during the secondary flare. No evidence 
for line emission was found with ~ 60 eV up- 
per limit on an iron line equivalent width. 

By combining the RXTE data with the 
Tcnma data, in't Zand, Baykal, and Strohmayer 
(1998) obtained a slightly improved orbital fit. 
This allowed a good determination of the or- 
bital period (8.3753 ±0.0003 days). However, 
the eccentricity is still not well determined, 
with a 68% confidence region of 0.14-0.38. 
They also reported the discovery of 18 s oscil- 
lations during a 1 hr flare, which may be the 
result of the beat frequency between the spin 



period and an accretion disk. In addition, 
the measured spin period Pg = 440.34llQQj^y 
suggests the pulsar has been steadily spinning 
down since 1983, with Pg = -|-0.225s yr^^. 

In this paper we present timing and spec- 
tral analysis of four ASCA observations of 
4U1907+09 and spectral analysis of an addi- 
tional 9 RXTE observations, roughly equally 
spaced around the 8.4 day orbit. We also re- 
port on archival observations of this source 
with HEAO A-1 and the RXTE All Sky Mon- 
itor. We then compare these and all previous 
observations of 4U1907-I-09 to a model of ac- 
cretion from a spherical wind coming from the 
companion. 

2. Archival Observations 

We have investigated archival HEAO A- 
1 scanning data on 4U1907-F09. HEAO A- 
1 observed 4U1907+09 in a scries of 32 s 
scans during 1978 April 7 - 14. The data 
show orbital behavior consistent with the 
Tenma/EXOSAT period with two peaks at 
orbital phases consistent with what was seen 
by Ariel V. The lightcurve from the A-1 ob- 
servations is shown in Figure 1. 

The RXTE All Sky Monitor has been ob- 
serving this source since it began regular ob- 
servations in February, 1996. We report here 
on ~ 1400 days of data (MJD 50136-51549). 
Errors on individual flux measurements are 
relatively large, making examinations on time 
scales of days unreliable. However, epoch fold- 
ing of the light curve clearly shows the orbital 
modulation. In Figure 2 we show the aver- 
age orbital light curve for the entire data set, 
with the orbital period and epoch from 
in't Zand, Baykal, and Strohmayer (1998). 
The vertical lines represent the 68% confi- 
dence region of the time of periastron. There 
is evidence for a secondary flare at a phase of 
~ 0.15. We also folded separately the light 
curves from the 3 energy channels (1.3-3.0 
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Fig. 1. — Light curve of 4U1907+09 during 
the 1978 HEAO A-1 observations. The verti- 
cal Unes in this and the following figures rep- 
resent the 68% confidence region of the time 
of periastron. 

keV, 3.0-5.0 keV, and 5.0-12.1 keV) of the 
ASM and found evidence for the secondary 
flare at a consistent phase in each. A com- 
parison of the first half of the data set to the 
second showed this is a fairly stable feature. 
The ratio of the flare to quiescent flux is dif- 
ficult to estimate since there is known to be a 
positive flux bias in the ASM data. In the case 
of 4U1907-I-09 this bias is exacerbated by sev- 
eral nearby sources including GRS 1915-1-105 
and W49B. 

We have also performed a period search on 
the ASM data in order to search for other pe- 
riodicities such as the proposed 41 day peri- 
odicity (Priedhorsky and Terrell 1984). We 
binned the data into 20 phase bins and cal- 
culated Xu for a constant model averaging 
over phase, as recommended by Collura et 
al. (1989). As Figure 3 shows, the 8.375 
day orbital period and several multiples of 
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Fig. 2.— RXTE All Sky Monitor light curve 
of the 8.3753 day orbital period averaged over 
the 4 year data set. The dotted line is the level 
of positive flux bias assumed for later model 
fits. 

it are poorly fit by a constant model. Al- 
though not evident in the ASM data fold- 
ing, five times the orbital period would cor- 
respond to 41.9 days which is consistent with 
the quasi-periodicity seen in the Vela 5B data 
by Priedhorsky and Terrell (1984) and inter- 
preted as an additional, longer period. We 
find no evidence for a periodicity in the long- 
term light-curve of this system which is unre- 
lated to the orbital period. 

3. ASCA Observations 

In August 1996, ASCA made four 10 ks 
observations of 4U1907-I-09, roughly equally 
spaced throughout the 8.4 d orbit. Observa- 
tional parameters are listed in Table 1. Ex- 
traction regions of 6' radius were used for the 
GIS detectors, while a 4' radius was used for 
SISO. A slightly smaller region was used for 
SISl (~ 3.5') since the source was nearer the 
edge of the chip. Since the source is in the 
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Fig. 3.— Period search of the RXTE Ah Sky 
Monitor data showing the 8.4 day orbital pe- 
riod and multiples. 

galactic ridge a local estimate of the back- 
ground is necessary. Therefore, source free 
regions of the detectors were used for extract- 
ing background light curves and spectra with 
a background rate of ~ 3 x 10~^cts/s in each 
detector. 

3.1. Temporal Analysis 

The background subtracted GIS3 light 

curve for all 4 observations (Aobsl-4) is shown 
in Figure 4. Aobsl consists of short flaring 
episodes with times of low flux in between 
where the count rate dropped by a factor of 
~ 10, reminiscent of the "dipping" seen in 
the earlier RXTE observations (in't Zand, 
Strohmayer and Baykal 1997). Aobs2 and 
Aobs3 show the source in a very low state, a 
factor of 50 below the average flux in Aobs4 
and about a factor of ~ 4 below the low flux 
states in Aobsl. Since Aobs4 has the best 
statistics and is relatively constant on time 
scales longer than the pulse period, we used it 
to find the best pulse period of 440.35 ± 0.06 



sec (after correcting the arrival times for the 
orbital motion of the pulsar, assuming the 
orbital parameters of in't Zand, Baykal, and 
Strohmayer (1998)) by epoch folding the light 
curve. We then folded each of the observa- 
tions at this period and show these in Fig- 
ure 5. The pulses seem to change shape and 
become weaker during the low states (Aobs2 
and Aobs3). 
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Fig. 4.— The ASCA GIS light curve of the 
four observations of 4111907+09 (upper panel) 
and of just Aobsl (lower panel). 

3.2. Orbital Phase Spectroscopy 

Using Xselect, spectra were extracted from 
each of the four detectors for each observa- 
tion. These were analyzed using XSPEC and 
fit to an absorbed power law model, adding 
narrow (a = 0) Gaussian lines if they signif- 
icantly improved the fit. Aobs2 and Aobs3 
had very low statistics so the two GIS de- 
tectors were combined into one data set, as 
were the two from the SIS detectors. In or- 
der to use all available photons, SIS Bright 
mode data were used. The results are shown 
in Table 2. It can immediately be seen that 
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Fig. 5. — Pulse profiles from the ASCA ob- 
servations. The photon arrival times are cor- 
rected for the pulsar orbit, and folded from 
a common epoch using the best fit period of 
440.35s from epoch folding obs4. The error in 
the period results in a relative phase uncer- 
tainty between obsl and obs4 of ~ 0.15. 

the lower states (Aobs2 and Aobs3) show no 
evidence of increased absorption. The lowest 
absorption occurs during Aobs4, the brightest 
observation, and is a factor of 2 lower than the 
high absorption seen during Aobsl. The spec- 
tral slope changes by a small, but statistically 
significant, amount. Weak iron line emission 
is present during Aobs4 and may be present 
during the other observations, although the 
statistics were not good enough to give more 
than upper limits at the 90% confidence level. 

In order to test whether the drops in flux 
during Aobsl are due to an increase of col- 
umn depth, we extracted spectra separately 
from the flaring and the low portions of the 
data set. As shown in Figure 6, there is no 
evidence of increased absorption during dips. 
In fact, there is evidence that the absorption 



column decreases during the dips, although 
the statistics on the low flux are such that 
the data are barely consistent within the 90% 
multi-parameter confidence region. There is 
no evidence of a significant change in photon 
index. If we assume the photon index remains 
constant, the change in absorption column is 
~ 2x lO^^cm"^. In the SIS data during the pe- 
riods of low flux there appears to be an emis- 
sion feature near 3.8 keV, which is evident in 
both SISO and SISl, although not apparent in 
the GIS data. In the GIS data the fit could 
be significantly improved by an inclusion of a 
width Gaussian line at 2.2 kcV. This feature 
appears in several of the data sets, although 
since it is near the known bump in the re- 
sponse due to Gold in the telescopes, this fea- 
ture should be viewed with caution. 
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Fig. 6. — High state and low state SIS energy 
spectra for Aobsl. 

During Aobs4 the source was only varying 
by a factor of 2 or so, mostly due to the pulsa- 
tions, and had a fairly high count rate allow- 
ing tight constraints on Nh and the photon 
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index. The fits are significantly improved by 
the inchision of an iron hnc at ~ 6.4 keV. The 
count rates were sufficiently high during this 
observation to extract spectra as a function 
of pulse phase using the best fit pulse period. 
There is no strong evidence for variability in 
Nfj- A strong anti-correlation between pho- 
ton index and intensity is evident which we 
show in Figure 7. There also seems to be 
a concentration of iron line emission in the 
pulse phase 0-0.1 spectrum. However, the un- 
certainties in the fit line strength at the other 
phases preclude a definitive measurement of 
line strength variability. 
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Fig. 7. — The photon spectral index as a func- 
tion of pulse phase (arbitrary epoch) during 
Aobs4. Errors represent the 90% confidence 
region. The dotted histogram is the GIS count 
rate. 



4. RXTE Observations 

RXTE observed 4U1907+09 for a second 
series of observations during December 1996. 
A total of 9 observations (Robsl-9), each ap- 
proximately 2 ksec in duration, were made 
spaced roughly one day apart in order to 



evenly sample the 8.4 day orbit. Data from 
the Standard2 observational mode were an- 
alyzed using FTOOLS 4.2. Combined spec- 
tra were extracted from all 5 PCUs (except 
obsl, which only had PCUs 0-3 on during 
the observation) and all layers. Model back- 
ground spectra were then constructed us- 
ing the background models provided by the 
RXTE Science Operations Facility, which take 
into account the particle induced background. 
However, they do not take into account the 
X-ray background from the galactic plane. 
To construct an accurate background model 
for 4U1907+09, it is necessary to take into 
account contributions from both the Galac- 
tic ridge and the nearby supernova remnant 
W49B. During two of the observations, Robs4 
and Robs5, the source was in the "dip" state. 
The model background was subtracted from 
these observations, and their spectra were 
examined. The count rate of ~ 18cts/s is 
comparable to what is observed by RXTE 
for the diffuse emission from the Galactic 
ridge (Valinia and Marshall 1998). We cre- 
ated a crude model of the Galactic ridge plus 
SNR W49B emission consisting of two ther- 
mal bremsstrahlung sources with fixed tem- 
peratures of l.SkeV and ~ 10 keV, moderate 
(1.6 X lO^^cm"^) absorption, and a ~ 6.7 keV 
gaussian lines with only the normalizations 
allowed to vary. These parameter values were 
taken from Koyama (1989) for this position in 
the galactic ridge and Smith et al. (1985) for 
W49B. This model adequately fit the data for 
Robs4 and Robs5 with normalizations that 
could reasonably be expected if the source 
had dropped to flux levels consistent with the 
low states seen in the ASCA data and the 
spectra were dominated by emission from the 
Galactic ridge and the nearby thermal super- 
nova remnant. Therefore, these observations 
were used to make a background spectrum for 
use in spectral fits of the other observations 
in the following way. Using the background 
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models available from the RXTE GOF, we 
generated background spectra for all 9 ob- 
servations. These were directly subtracted 
from the observation spectra using the math- 
pha FTOOL to create the files to be used in 
the fitting. The background subtracted pha 
files of Robs4 and Robs5 were then added to- 
gether to create a local X-ray background for 
4U 190 7+ 09 which was used as the background 
file in XSPEC. 

All the spectra from the remaining seven 
observations were modelled using an ab- 
sorbed power-law with a high-energy cutoff 
in XSPEC. The data were fit in the energy 
range of 2.5 - 20 keV, and the results are 
listed in Table 4. The measured cut-offs are 
consistent with the broad-band spectral fits 
of Beppo-Sax (Cusumano et al. 1998). We 
also obtained the archival RXTE data from 
February 1996 first analyzed by in't Zand, 
Strohmayer and Baykal (1997) to fit for av- 
erage fluxes and absorption in the 2.5-10 keV 
band for inclusion in the modelling below. 

5. Modelling of Orbital Spectral Vari- 
ations 

In Figure 8 we show measurements of nH 
as a function of orbital phase for spectral 
measurements from ASCA and RXTE (this 
work), EXOSAT (Cook and Page 1987), 
Tenma (Makishima et al. 1984), and SAX 
(Cusumano et al. 1998). As a consistency 
check, we also plot the ASM hardness ratio 
of (5.0-12.1 keV cts/s)/(1.3-3.0 keV cts/s) as 
a function of phase (dotted line). Since the 
spectral slope does not vary greatly, we ex- 
pect this hardness ratio to change mostly due 
to column depth. The vertical lines represent 
the 68% confidence region of the time of peri- 
astron. We see that the absorption begins to 
increase towards the end of the primary flare 
stays high until right before the secondary 
flare, where it drops again. In Figure 9 we 
plot the photon spectral index as a function 
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Fig. 8. — The multi-mission measurement 
of nH as a function of orbital phase. Mea- 
surements are from ASCA (open triangles), 
Feb. 1996 RXTE (filled triangles), SAX (open 
square), EXOSAT (filled square), Tenma 
(open hex), and Dec. 1996 RXTE (filled hex). 
The histogram is the ASM hardness ratio of 
the countrates from channel 3 over channel 1. 

of orbital phase. Small but significant changes 
are observed. In Figure 10 we plot the 2-10 
keV flux measurements from these same ob- 
servations versus orbital phase with the ASM 
lightcurve overplotted. The secondary flare 
is twice seen to be nearly as bright as the 
primary during these observations. One of 
the RXTE observations seen at an overlap- 
ping phase with the ASCA Aobs2 is still seen 
to be moderately bright [Fx ^ 10^^'^ ergs 
cm~^ s~^) although dipping activity is seen 
during the observation. Therefore, extended 
periods of low flux as seen in the ASCA Aobs2 
and Aobs3 probably do not occur every or- 
bit. In general there is a positive flux bias 
of ^ 0.1 cts/s in ASM data. However, in 
crowded regions of the Galactic plane, such 
as around 4U1907+09, there arc likely to be 
further systematic biases. By comparing the 
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folded ASM light curve to the observed multi- 
mission fluxes we estimate a positive flux bias 
of ~ 0.2 — 0.6 cts/s. We will assume a bias of 
0.3 for the analysis below. 
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Fig. 9. — The multi-mission measurement 
of the photon index as a function of orbital 
phase. 

In this section we model the primary flare 
and absorption as accretion from a spherical 
stellar wind with the addition of an absorb- 
ing gas stream. This model, first applied to 
Tenma observations of GX 301-2, is detailed 
in Leahy (1991). We first attempt to model 
the average ASM lightcurve and multi-mission 
absorption data as accretion from a simple 
spherical stellar wind with velocity profile: 



'^wind 



.(1 - Rc/r) 







where Voo is the terminal velocity, assumed to 
be 3 times the escape velocity, and Rc is the 
radius of the companion. The Castor, Abbott, 
and Klein (1975) theoretical wind model has 
P = 0.5, but observations of OB stars (Lei- 
therer 1988; Boroson ct al. 1999) suggest 
that P ~ 0.7 — 2 depending on spectral type. 
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Fig. 10. — The multi-mission measurement 
of the 2-10 keV flux cis db function of orbital 
phase. The histogram is the folded lightcurve 
from the ASM, with the y axis starting at the 
assumed flux bias of 0.3 cts/sec. 

We assume (3 = 1, consistent with the as- 
sumption of an OB supergiant. In the case 
of a lower mass Be-type star, models of the 
wind structure are more complex and uncer- 
tain. However, the velocity profile of the polar 
wind from a Be-star is similar to that of an OB 
supergiant (de Araujo, de Freitas Pacheco, 
and Petrini 1994), and the profile of the equa- 
torial wind probably has a similar functional 
form but lower terminal velocity. If the com- 
panion is a Be-star, it is likely that the spin 
axis is not aligned with the orbital axis (see 
discussion) and so emission dominated by ac- 
cretion from the polar wind over much of the 
orbit may not be an unreasonable assumption. 

We assume the orbital period P = 8.3753d, 
longitude of periastron (i.e. the observer's 
viewing angle) to = 330 dcg, and periastron 
phase of the best fit model from in't Zand, 
Baykal, and Strohmayer (1998) and fit the 
ASM light curve for a range of inclination an- 
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gles by adjusting the eccentricity and the pa- 
rameter X = M2o/'^3) where M20 is the mass 
loss rate in units of lO^'^ grams/sec and is 
the distance in units of 3 kpc. We choose as 
our canonical values i = 47.9, which implies a 
companion mass Mc = 24Mq given the mass 
function derived from the orbital fits assuming 
a neutron star mass of 1.4M0, and i = 75.7° 
for a companion mass Mc = 12Mq more in 
line with a Bc-star. 

For the companion radius we follow van 
Kerkwijk, van Oijen, and van den Heuvel 
(1989) and assume Rc = 31Rq for a 24Mq 
star. However, a star of this radius would be 
larger than its Roche lobe at pcriastron, as 
defined in Brown and Boyle (1984), assum- 
ing the best fit eccentricity of e = 0.28 (in't 
Zand, Baykal, and Strohmaycr 1998), a ro- 
tational velocity Vsini = 90 kni/s (van Kerk- 
wijk, van Oijcn, and van den Heuvel 1989), 
and alignment of the rotational and orbital 
axes. In Figure 11 we plot the Roche lobe ra- 
dius as a function of companion mass as well 
as the limits on radius inferred from the lack 
of eclipses. All supergiant companions would 
greatly overflow their Roche lobes and would 
have to be excluded if we assume this 
dius limit. However, a non-rotating (or non- 
aligned) star has a greater Roche lobe radius, 
as would a system where the eccentricity is 
less. For the 12Mq star, we use the lyc (1986) 
value of Rc = 6Rq for a Be-star compan- 
ion. For intermediate inclinations and masses, 
the apropriate radius is dependent on spectral 
type. For supergiants, the radius actually in- 
creases with decreased mass, which would ex- 
acerbate the Roche lobe problem. However, 
a type lb, HI, or V are all smaller than our 
canonical la supergiant star. In Figure 11 we 
plot radii for a number of early B stars (Un- 
derbill and Doazan 1982). A linearly interpo- 
lation between our two canonical stars seems 
to provide a reasonable mass-radius relation- 
ship for modelling purposes. 
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Fig. 11. — Assumed Mass/Radius relation- 
ship (solid line) compared to eclipse limit and 
periastron "Roche Lobe" radii assuming dif- 
ferent eccentricities and stellar rotation rates 
in the orbital plane (O*). Points are various 
stellar values for main sequence, giant, and su- 
pergiant B0-B3 stars (Underhill and Doazan 
1982) 

In Figure 12 we show the range of model 
curves compared to the ASM light curve. To 
fit the higher mass companion to the shape 
of the primary flare an eccentricity e < 0.14 
is required (which would also allow the star 
to underfill its Roche lobe at periastron). For 
the lower mass companion the best fit to the 
observed primary flare peak is e = 0.26, more 
in line with the best fit eccentricity from the 
pulse arrival analysis. In no case can a sec- 
ondary flare be produced from this model. 

This luminosity model implies a varying 
absorption function of orbital phase, 

with the overall normalization being only a 
function of M. Using the parameters given 
above and an assumed interstellar absorp- 
tion of nH = 1.5 X lO^^cm^^, we find that 
-^20 = 1 — 1-5 best fit the absorption data over 
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Fig. 12. — Spherical wind model fits to ASM 
lightcurve assuming inclination angles of 47.9° 
(Mc = 24M0), 60° (Mc = I6M0), and 75.7° 

(Me = I2M0). 

our range of inclination angles. This is some- 
what low for a B supcrgiant, but quite high for 
a Be-star (de Araiijo, de Freitas Pacheco, and 
Petrini 1994). In order to improve the fit, 
wc also vary the the longitTidc of periastron 
since that is poorly constrained by the pulse 
arrival time analysis of in't Zand, Baykal, and 
Strohmayer (1998). In Figure 13 we plot the 
results compared to the multi-mission data. 
Neither the phase or the amplitude of the 
large increase in absorption after periastron 
can be produced from this simple spherical 
wind model. 

We can use the M from the absorption 
fits and the parameter x from the luminos- 
ity model to estimate distances to the system. 
For the 2AMq model, we derive a distance 
d ~ 4.3 kpc, while for the WMq and I2M0 
models, we derive d ~ 1.0 kpc and d ~ 0.4 
kpc respectively. The first is in line with the 
distance estimates of 5.9 to 2.4 kpc for a su- 
pergiant companion from van Kerkwijk, van 
Oijen, and van den Heuvel (1989), and the 
latter are consistent with their estimates of 
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Fig. 13. — Spherical wind model absorption 
compared with multi-mission spectral fits. 

1.5 to 0.6 kpc for a giant. Note that the min- 
imum measured nH ~ 1.5 x lO^^cm"^ means 
the assumed interstellar component of the ab- 
sorption is an upper limit. Thus the derived 
values of M, and hence the inferred distances, 
should be considered lower limits. 

To fit the post-periastron absorption, we 
added a linear trailing gas stream to the above 
model (see Leahy (1991) for details). Such a 
stream can be interpreted as arising from the 
wake caused by the passage of the neutron 
star through the wind, or from an accretion 
stream coming from the companion (although 
this latter interpretation would require modi- 
fication of the luminosity model). The stream 
is modelled as a density enhancement with a 
gaussian profile perpendicular to the stream 
line, an exponential cutoff on the front side 
and no cutoff in the downstream direction. 
The total contribution to the absorption is 
then obtained by integrating along the line of 
site. For a given orbital phase, longitude of 
periastron, and orbital inclination , the am- 
plitude of the absorption enhancement of the 
trailing gas stream is completely determined 
by the prodiict UoCJ ^ lO^^cm"^, where Ho is 
the central density and a the gaussian width. 
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The additional absorption due to the stream 
is shown in the top of Figure 14. The phase of 
the absorption peak is somewhat earher than 
indicated by the data, so we also tried offset- 
ting the stream line towards the companion 
star by angles of 20° and 40°. The bottom 
panel of Figure 14 shows how the combined 
wind plus stream model compares to the ab- 
sorption data for the various cases. We find 
that an adequate fit can be obtained with the 
larger inclination angles and an offset stream. 
A downstream cutoff of the stream would tend 
to broaden the peak somewhat, which might 
further improve the fit. 




0.2 0.4 0.6 0.8 1 



Orbital Phase 

Fig. 14. — Absorption due to a trailing gas 
stream. The top panel shows the orbital vari- 
ation of absorption due to the stream for vari- 
ous parameters. The bottom panel shows the 
wind plus stream absorption model compared 
to the multi-mission measurements of nH. 

6. Discussion 

If the eccentricity of e > 0.2 of the best 
fit orbital parameters from pulse arrival time 
analyses is correct, the primary flare seen in 
the overall light curve can be fit using a model 
of accretion from a spherical wind off of a star 
of mass Mc < 16Mq. However, the absorp- 



tion measurements and the persistence of the 
secondary flare require modifications to this 
model. We have shown that a trailing ac- 
cretion stream moderately offset towards the 
companion can account for the extra absorp- 
tion component, although not the secondary 
flare. Such an offset stream may be indicated 
by the high resolution two-dimensional hydro- 
dynamical simulations of Benensohn, Lamb, 
and Taam (1991). The "flip-flop" instability 
seen in their simulations modifies the direc- 
tion of the downstream flow so that it switches 
back and forth resulting in a wake that could 
be characterized by two offset streams, one of 
which could provide the enhanced absorption 
seen in the data (see Fig. 2 of Benensohn, 
Lamb, and Taam (1991)). This instability 
can in principle also explain the "dipping" or 
"flaring" behavior seen in Aobsl and in the 
RXTE observations of in't Zand, Strohmayer 
and Baykal (1997), and may even be compat- 
ible with a drop in absorption during low ac- 
cretion phases. However, this instability may 
just be an artifact of two-dimensional simu- 
lations, and whether similar affects occur in 
a three-dimensional flow is an open question. 
Recent 3D simulations by Ruffert (1999) do 
exhibit active unstable phases, although the 
changes in accretion rate tend to be smaller 
and the downstream density profile is unclear. 

The spherical wind accretion models favor 
a high inclination angle, implying a compan- 
ion mass more in line with a Be-star than 
a supergiant. In addition, if the companion 
mass is high, the models imply an eccentric- 
ity near the low end of the acceptable range. 
This lower eccentricity is also required to keep 
a massive supergiant from greatly overflow- 
ing its Roche lobe near periastron. Even with 
a smaller eccentricity, it is likely that peri- 
odic Roche lobe overflow (Brown and Boyle 
1984) or tidal stripping (Layton et al. 1997) 
would be the dominant mass transfer mecha- 
nism. In this case, the observed accretion rate 
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and orbital variation may be somewhat low, 
unless long-term storage of material by an ac- 
cretion disk could keep the accretion rate low 
except for occasional large outbursts like the 
one seen by Ariel V in 1980 (Marshall and 
Rickets 1980). If the best fit eccentricity 
value of e ~ 0.28 is correct, it would seem 
to exclude the accretion from a supergiant 
wind model. This suggests that an alternative 
model of the accretion perhaps should be con- 
sidered. One persistent possibility is an equa- 
torially enhanced wind with the neutron star 
in an inclined orbit, as suggested by Iye(1986). 
In this case, the neutron star passes behind an 
extended matter disk around the companion 
after periastron and comes back through af- 
ter apastron. While tidal forces would tend to 
align the neutron star orbital axis with that 
of the companion's spin axis, the timescale for 
this is similar to the circularization timescale, 
and accretion processes might coTintcract the 
tidal forces somewhat (Matese and Whitmirc 
1983). Given that the orbit is eccentric, we 
do not view a somewhat inclined orbital axis 
as being unreasonable. 

In the equatorially enhanced wind scenario, 
the flares are fed by matter gathered by the 
neutron star as it passes through the matter 
envelope. If the outflow speed is somewhat 
slower near the equator, the size of this en- 
velope might be limited by the neutron star's 
orbit. Such truncation is indicated from stud- 
ies of the Ha emission line in several Be/X-ray 
binary systems (Zamanov et al. 2000). The 
matter is depleted by the periastron passage 
of the neutron star, so significantly less ma- 
terial is available further out at the apastron 
radius. If sufficiently depleted, there may be 
no secondary flare at all. Hence we get an in- 
termittent but phase locked secondary flare. 

The major difficulty to this picture is the 
"virtually certain" classification of the com- 
panion as a supergiant by van Kerkwijk, van 
Oijen, and van den Heuvel (1989). This is 



based on a collection of observed properties of 
the system, with some favoring a Be-star, and 
some favoring a supergiant. Two factors seem 
to greatly favor the supergiant interpretation. 
First is the color excess Eb-v = 3.4 ± 0.2, 
which is quite high for the inferred distance of 
< 1.5 kpc of a Be-star. However, the variation 
in absorption column with orbital phase sug- 
gests that much of this could be local. Indeed, 
the average excess determined from interstel- 
lar absorption lines is only Eb-v = 2.4 ± 0.5 
or Eb~v = 2.7 ± 0.6, depending on which co- 
efficients for the ratio of equivalent width to 
color excess used (van Kerkwijk, van Oijen, 
and van den Heuvel 1989). These lower val- 
ues are more in line with the inferred interstel- 
lar X-ray absorption of nH < 1.5 x lO^^cm"^. 
The independent distance estimate van Kerk- 
wijk, van Oijen, and van den Heuvel (1989) 
obtain from the interstellar sodium lines of 
d = 2.1 lb 0.6 is between the supergiant and 
Be-star distances, and consistent with both 
given the error. 

The second factor strongly favoring the su- 
pergiant interpretation is the large FWHM~700 
km/s of the Ha line combined with the 
relatively low projected rotational velocity 
Vsini ~ 90kms. While neither value is ex- 
ceptionally unusual for a Be-star, the combi- 
nation is (Dachs et al. 1986). The inclination 
angle suggested by the observed projected ro- 
tational velocity is i ~ 15° assuming typical 
Be-star rotational velocities, as opposed to 
the inferred orbital inclination i ~ 60° — 75° 
of the neutron star with a giant or main se- 
quence star companion, consistent with the 
inclined orbit scenario. However, the large 
value of the FWHM suggests the outflowing 
disk is seen nearly edge on. The circumstellar 
envelope can and would be affected by the 
neutron star, as evidenced by the large-scale 
perturbations observed in the Ha profiles of 
several Be/X-ray binaries (Negueruela et al. 
1998). Since the orbital period of 4U1907-h09 
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is shorter than any of the known Be/X-ray 
binaries (except A0538-66, which is not weh 
observed), it would be expected to have a 
much greater perturbing effect. It is inter- 
esting to note that the projected orbital ve- 
locity of a particle in the orbital plane near 
the Roche lobe radius of a lower mass star in 
4U 1907+09 is roughly half the observed Ha 
FWHM. If the gravitational attraction of the 
neutron star could warp the equatorial out- 
flow of the companion such that it would be 
nearly co-planar to the neutron star orbit at 
the Roche lobe radius, the observed FWHM 
might result. Hydrodynamical simulations 
should be able to determine if such a scenario 
is feasible. 

The relatively long periods of low accre- 
tion rate seen in Aobs2 and Aobs3 are some- 
what puzzling. Inhomogeneities in the wind 
from the companion could possibly account 
for these variations, although the apparent 
duration of several days during the ASCA ob- 
servations would argue against this. If the 
accretion is modified by a transient accretion 
disk, instabilities in the accretion disk could 
account for both the sudden changes in ac- 
cretion rate seen in Aobsl and the longer pe- 
riods of quiescence seen in Aobs2 and Aobs3. 
The consistent spin-up of the neutron star and 
the 18 sec. oscillations seen by RXTE (in't 
Zand, Strohmayer and Baykal 1997) argue 
for the existence of an accretion disk. The on- 
off behavior seen by both RXTE and ASCA 
are suggestive of magnetically inhibited accre- 
tion from a disk whose inner edge is near the 
corotation radius. The residual emission dur- 
ing the off periods may then be a result of 
accretion on to the magnetosphere or a bow 
shock. However, the slow spin period and 
correspondingly large corotation radius makes 
this unlikely. 



7. Comparison with GX301-2 

4U1907-I-09 is similar in many respects to 
the wind fed Supergiant/X-ray pulsar binary 
system GX301-2. Both have eccentric orbits 
and occupy a similar location on the spin- 
period versus orbital-period diagram (Cor- 
bet 1986) (for GX 301-2, e ~ 0.47, Pspin 
680s, Por6 ~ 41. 5d, Koh et al. (1997)). 
The absorption profile and primary orbital 
flux peak of both systems have been success- 
fully fit by a spherical wind plus gas stream 
model (Leahy 1991). GX 301-2 also has ex- 
tended periods where there is a secondary or- 
bital flare at a consistent phase that cannot 
be fit by wind models. An equatorialy en- 
hanced outflow has also been suggested for 
this system and models of such an outflow can 
grossly fit the orbital light curve (Koh et al. 
1997). Both systems have companions with 
similar temperatures (spectral type B2) and 
projected rotational velocities (vsini ~ 70 for 
GX 301-2, Parkes et al. (1980)). 

However, the companion star of GX 301- 
2, Wray 977, is unambiguously identified as 
a supergiant by the strength of various lines 
in the blue portion of the spectrum (Parkes et 
al. 1980). The observations of lye (1986) and 
van Kerkwijk, van Oijen, and van den Heuvel 
(1989) did not extend into the blue portion 
of the spectrum, so such diagnostics are not 
available for 4U 1907-^09. The H/3 emission 
line of Wray 977 has a P-Cygni profile (Parkes 
et al. 1980), clear evidence of an expanding 
atmosphere. No such profiles have been seen 
in 4U 1907+09. 

The local absorption column and iron 
line equivalent width of the X-ray spectrum 
of GX 301-2 are both more than an order 
of magnitude stronger than in 4U1907+09. 
This implies the mass loss from Wray 977 
is much greater than from the companion of 
4U1907+09, despite the similarity in the pro- 
posed spectral classes. While enhanced mass 
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loss in these systems could be expected due to 
gravitational interaction with a neutron star, 
the orbit of GX 301-2 is significantly wider, 
and there is no question that the companion 
remains inside its Roche lobe (although at pe- 
riastron it can come close, Leahy (1991)). If 
both were supergiants, the naive expectation 
would be greater mass loss from the com- 
panion of 4111907+09 due to transient Roche 
lobe overflow (Brown and Boyle 1984) or tidal 
stripping Layton et al. (1997). 

The changes in spectral index with pulse 
phase of both systems show a softening of the 
spectrum near the pulse minimum in the 1-10 
keV band (Saraswat et al. 1996; Leahy and 
Matsuoka 1990). The iron line equivalent 
width in GX 301-2 varies with pulse phase as 
well, with a peak near the pulse minimum. 
The possible iron line enhancement near the 
pulse minimum in the 4U1907+09 ASCA 
data, if real, would be consistent. These sim- 
ilarities in the phase-resolved spectra suggest 
the accretion columns are similar. This is 
not surprising, since the spin periods and sur- 
face magnetic field strengths (as inferred from 
their cyclotron absorption lines) are similar. 

8. Conclusions 

Observations of the variation of the spectra 
of 4U1907+07 with orbital phase are consis- 
tent with a simple spherical wind accreting 
onto the neutron star from an OB supergiant 
companion with excess absorption caused by 
a trailing accretion wake. However, the im- 
plied inclination angle for such a model sug- 
gests a relatively low mass (< 16Mq) for the 
companion. The excess of absorption column 
seen after periastron might also be interpreted 
as the neutron star passing through an equa- 
torially enhanced matter envelope in an in- 
clined orbit. Observations of the secondary 
fiare near apastron by RXTE and HEAO A- 
1, as well as the average orbital light curve 
from the RXTE ASM, further support the hy- 



pothesis that this flaring has been occurring 
at a consistent phase but with variable ampli- 
tude for 20 years. The measurements of the 
low states by ASCA suggest that the accre- 
tion rate can change by 1-2 orders of magni- 
tude over a short time period, perhaps due to 
a transient accretion disk. The most critical 
components to a better understanding of this 
system are more secure measurements of the 
orbital eccentricity and periastron phase and 
improved optical measurements of the com- 
panion. A high resolution optical spectrum 
at shorter wavelengths and orbital radial ve- 
locity measurements could unambiguously de- 
termine the companion type. 
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sions as to the nature of the companion star 
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Table 1 

ASCA Observation Parameters 



Observation 


Sequence ID 


Date 


Orbital Phase 


On-Time 


Mean CIS rate 






MJD-50300 


= 


ks 


Cts/s 


Aobs 1 


44020040 


68.65-69.06 


0.93-0.98 


10.2 


0.702 ± 0.008 


Aobs 2 


44020010 


70.75-71.12 


0.18-0.22 


10.1 


0.067 ± 0.005 


Aobs 3 


44020020 


72.34-72.53 


0.37-0.39 


9.3 


0.064 ± 0.005 


Aobs 4 


44020030 


74.54-74.94 


0.63-0.68 


9.9 


3.447 ± 0.015 



Table 2 
ASCA Spectral Fits 



Observation 




r 


line E 


line EqW 


F2-lQkeV 


X^(dof) 




xlO^^cm"^ 




keV 


eV 


xl0~^° ergs cm~^s~^ 










Fits from GIS2 


+ GIS3 






Aobs 1 


3.89l°;i 


1 02+°-l2 

J--U^-0.12 


fi 53+°-32 


76 ±33 


1.11 ±0.02 


0.58(38) 


Aobsl hi 


4.181^:1^ 




6.45l°:i 


82 ±33 


1.73 ±0.02 


0.50(38) 


Aobsl lo 


r, ir, + 1.68 

^•-'■^-1.12 


0.92;°;i 






0.23 ± 0.02 


0.24(40) 


Aobs 2 


2.5911:^^7 


2 07+°-^^ 

-0.65 






0.066 ± 0.005 


0.25(40) 


Aobs 3 


2.70l}i^ 


1 nQ+0.88 
-^•^•3-0.73 






0.066 ± 0.005 


0.19(40) 


Aobs 4 


1 74-^0.06 

-•-'^-0.06 


-, 90+0.04 
^•^■^-0.03 


6.53l°- 
Fits from SISO 


79± 16 
± SISl 


4.20 ± 0.02 


1.71(38) 


Aobs 1 


q Oq+0.39 

"^•O'J-0.36 




6.34 (frozen) 


< 57.1 


0.87 ± 0.01 


0.53(289) 


Aobsl hi 


4ni+0.42 




6.34 (frozen) 


< 65.8 


1.42 ± 0.02 


0.51(289) 


Aobsl lo 


-1.20 








0.18 ±0.01 


0.35(35) 


Aobs 2 


^•^"-0.93 


1 8Q+0.78 






0.051 ± 0.003 


0.30(28) 


Aobs 3 


1 oc;+1.23 
-'-•°^-0.84 


1 Q4+0.72 






0.063 ± 0.004 


0.19(28) 


Aobs 4 


-■-"'-o.oe 


1 qcr-l-0.04 


fi Q4+0.08 


70.0 ± 15.8 


3.49 ± 0.02 


0.83(324) 
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Table 3 

RXTE Observation Parameters 



Observation 


Obs ID 


Date 


Orbital Phase 


On-Time 


Mean raw PCA rate 






MJD-50400 


= 


ks 


Cts/s 


Robs 1 


10154-02-01 


36.40-36.45 


0.02 


2.4 


243 


Robs 2 


10154-02-02 


37.39-37.42 


0.14 


2.4 


238 


Robs 3 


10154-02-03 


38.39-38.42 


0.25 


2.4 


233 


Robs 4 


10154-02-04 


39.39-39.42 


0.37 


2.4 


181 


Robs 5 


10154-02-05 


40.39-40.42 


0.49 


2.4 


179 


Robs 6 


10154-02-06 


41.05-41.08 


0.57 


2.5 


291 


Robs 7 


10154-02-07 


42.06-42.12 


0.70 


2.6 


331 


Robs 8 


10154-02-08 


43.06-43.09 


0.81 


2.4 


333 


Robs 9 


10154-02-09 


44.05-44.09 


0.93 


3.3 


216 



Table 4 

RXTE Spectral Fits from Dec. 1996 Observations 



Observation 


Phase 


riH 
xlO^^cm"^ 


r 


Cutoff Energy 
keV 


F2-WkeV 

xlO~^° ergs cm~^s~^ 




Robsl 


0.995 


6 -ii+o-^s 

"•^-'^-0.43 


-, r,r 1+0.029 
±.Z01_g Q29 


1 q r9+0.25 
10.0Z_Q 25 


2.92 


0.63 


Robs2 
Robs3 


0.114 
0.233 


9 QQ+0.49 
^•<JO_o,49 


1 9c;c:+0.038 
-L-^<J'J_0.039 
1 9trq+0.067 


i4.ool°:S 

13.92lJ:^2 


1.63 
0.89 


0.84 
0.92 


Robs6 
Robs7 


0.550 
0.671 


2-68l°:i 

q 9Q+0.36 
^•^"-0.35 


1 IS/1+0-022 

1 Of5q+0.028 
1.0Dy_g Q29 


IS 47+0.21 

1 ^ t;9+0.28 
±O.OZ_Q 29 


3.30 
4.41 


1.00 
1.50 


Robs8 
Robs9 


0.790 
0.909 


7.89^°J^ 
7.9llgj^ 


-| 4-|t;+0.020 

1 QR9+0.049 
l.^DZ_Q Q4g 


13.50;gi? 
13.62l0i^ 


4.29 
1.21 


0.88 
0.72 
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